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ABSTRACT 

Recent Spitzer observations have revealed a substantial population of z ~ 2 
ULIRGs with deep silicate absorption (rg 7 >l). This paper reports a 20cm 
radio study of such a sample to elucidate their physical nature. We discover 
that a substantial fraction (40%) of deep silicate absorption ULIRGs at z ~ 2 
are moderately radio- loud with Li.4GHz — 10 25 - 10 26 WHz -1 . This is in strong 
contrast with z <1 radio galaxies and radio- loud quasars where none of the sources 
with available IRS spectra have r 9 7 > 1. In addition, we observe radio jets in two 
of our sources where one has a double lobe structure ~200kpc in extent, and 
another shows a one-sided jet extending ~ 90 kpc from the nucleus. The likely 
high inclination of the latter, coupled with its deep silicate absorption, implies 
the mid-IR obscuration does not share an axis with the radio jets. These sources 
are highly obscured quasars, observed in the transition stage after the birth of 
the radio source, but before feedback effects dispel the ISM and halt the black 
hole accretion and starburst activity. 

Subject headings: galaxies: active — galaxies: high-redshift — galaxies: jets - 
infrared: galaxies — radio continuum: galaxies 



Introduction 



The local relationship between the bulge mass and black hole mass (jMagorrian et al. 



19981 ) suggests that the stellar mass build-up and black hole growth were co-eval. This is 



supported by the num ber density evolution of quasars and radio galaxies (jWall et al.l 12005 



Richards et al.ll2006af) being similar to th e evolution of the star-formation rate (SFR) density 



with redshift (e.g. iGiavalisco et al.l 12004 ). Observations of significant far-IR e mission poten- 



tially suggest starbursts associated with high- z quasars, and radio galaxies (e.g. IWillott et al. 



20021 : iReuland et al.ll2003l : iBeelen et al.ll2006l ) and indicate rapid growth of the host galaxy 
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likely fueled by mergers (IHutchings et al.ll2006l ; iHopkins et al.l 120061 ). Studi es of the young 



stella r populations of radio galaxies also suggest merger-driven starbursts (ITadhunter et al. 



20051 ). Observations o 



the surrounding ISM (INesvadba et al.l 120061 ) 



z ~ 2 radio galaxies indicate the powerful outflows that will disperse 



To study the process of a quasar emerging from its dusty cocoon, we need transition 
Ultr aluminous Infrared Galaxies (ULIRGs) - i.e. containing a still dust enshrouded quasar 
(e.g. ICanalizo fc Stockton! 120011 ). Unlike infrared, optical, or X-ray diagnostics, extreme 
radi o luminosity is a dust-independent, unambiguous indicator of a highly obscured AGN 



(e.g. 



Martmez-Sansigre et al.ll2005l ; iDonley et al.l 120051 ) 



We have a sample of Spitzer-selected z ~ 2 UL IRGs whose con t inuum-dominated mid- 
IR spectra already imply the presence of an AGN (jYan et al. 2007 ; Sajina et al. 2007 ). In 
this letter, we highlight the seven sources of this sample which are radio-loud, confirming 
the AGN, while their deep 9.7/im silicate absorption indicates significant obscuration. 



Throughout this letter we adopt the ACDM cosmology with Qm 
J Ho = 70kms- 1 Mpc~ 1 . 



0.3, Q A = 0.7, and 



2. The sample 



Our sources are part of a sample of 52 high-z ULIRGs with IRS spectra (jYan et al. 



20071 ). In summary, the sample is selected from the Spitzer First Look Survey (FLS). They 
are selected for having 24 /im flux density brighter than 0.9mJy, and red 24-to-8 yum and 
24-to-R colorfl Of these, 47 have IRS-based redshifts, with the majority (33/47 = 72%) at 
1.5 < z <3.2. The analysis of the IRS spectra, and derivation of PA H luminositie s , equi valent 
widths, and silicate feature optical depths^] (T9.7) is presented in Sajina et al. ( 2007 ). All 
sources in our sample show evidence of both star-formation and AGN activity; however, based 
on their PAH-strengths vs. mid-IR continua, the 'weak-PAH' (EW7.7 < 0.8/im ) sources are 
believed to be energetically dominated by their AGN. About 60% of the these weak-PAH 
sources have significant 9.7 /im silicate absorption with 79.7 > 1. Roughly half of these show 
extreme optical depths (79.7 > 3). At such high obscurations the IR spectra of AGN an d 
compact, nuclear starbursts are indistinguishable (IRoussel et al.ll2006l ; iLevenson et al.ll2007l ). 



fl(24,8) = log 10 (^(24 M m)/^/ i/ (8/xm) > 0.5; (2) #(24,0.64) = Iog M (i// v (24/im)/i// tf (0.64/im) > 

1.0. 

2 Our 9.7 /mi optical depth is measured from a n extinction-correc ted continuum and relates to the observed 
depth of the silicate feature as r 9 7 = 1.4 rg, (see lSaiina et al.ll2007l for details). 
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Of the total sample, 60% are detected in the VLA 20 cm map of t he FLS (1 Condon et aU20031) . 
and 48% are detected in the GMRT 610MHz map of the field JGarn et al.ll2007h . For a full 
description of the radio data see Sajina et al. (2007, in prep.). 

In this letter, we only look at the 28 z>1.5 continuum-dominated (i.e. likely AGN 
dominated) sources. We use their radio emission to reveal the hidde n AGN. We take 
Liaghz > 10 25 WHz -1 as our 'radio-loud' criterion (e.g. IJiang et al.l 120071 ). while we discuss 
the more standard radio-to-optical ratio criterion in § 13.41 The seven sources which satisfy 
this criterion are li sted in Table [J along with a summ ary of their mid-IR properties from 
our earlier papers (jYan et al.l 120071 ; ISajina et al.l 120071 ) as well their radio fluxes and lumi- 
nosities from Sajina et al. (2007, in prep). These seven sources represent 25% of the z~2 
continuum-dominat ed sources. For c omparison, optically-selected quasars have ~10-20% 
radio- loud fraction (I Jiang et al.ll2007l ). 



3. Analysis 

3.1. Mid-IR obscuration 

The top panel of Figure [1] shows that all of our radio-loud sources have deep silicate 
absorption features (rg.7~l-6). Among our z~2 high-r sources, ~40% are radio-loud, 
making this an exceptionally efficient means of selecting high-z radio-loud sources. By con- 
trast, none of the low-r sources are radio-loud. The observed L 14Atm of the low- and high-r 
sources are not significantly different (this holds for Lfir as well; Sajina et al. 2007 in prep.). 
However, the extinction corrected 14 fim continua are significantly higher for the deeply ab- 
sorbed sources (bottom panel in Figured]). Higher mid-IR luminosities might indicate any 
of greater black hole masses, greater accretion rates, or more compact dust geometry lead- 
ing to hotter dust. At present it is unclear which, if any, of these reasons mig ht explain 



the la ck of radio-loud sources among our low-r continuum-dominated sources. IShi et al. 



(120061 ) show a correlation between the gas column densities (A%) and rg.7. Extrapolating 
this to the optical depths of our sources, we find Nn > 10 27 cm~ 2 . If confirmed, these sources 
qualify as Compton-thick quasars. It is however likely that this trend breaks down due to 
differences in the properties and geometry of the obscuring medium (for example a more 
uniform, screen-like obscuration as opposed to a more clumpy one). Initial Chandra obser- 
vations support strong X-ray obscuration (Bauer et al. 2008 in prep) given non-d etections in 



30 ks exposures. Lastly, five of the seven meet the lMartmez-Sansigre et al.l (120051 ) Compton- 
thick quasar criteria (MIPS15880 has too much 3.6 /im flux, while MIPS16122 has too much 
1.4GHz flux). 
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3.2. Spectral indices 



The radio spectral indices^] (a) based on the integrated fluxes at 610 MHz and 1.4GHz 
are listed in Tabled! The different resolution of the two images means that the 1.4 GHz 
map might be missing extended flux relative to the 610 MHz. We estimate the strength 
of this effect by comparing the a values based on the integrated or peak 610 MHz fluxes. 
We find that our a's can be overestimated by up to ~0.2dex. This mean up to 0.1 dex 
overestimation of the radio luminosities, which is not significant for our conclusions. The 
largest effect is found for our strongest radio source, MIPS16122, where a ranges between -0.0 
and -0.5 for the peak and integrated flux values respectively. This is our only potentially flat- 
spectrum source. The typ ically steep spectra we find are consistent with the a~-l of z~2 
Type-2 quasars found by iMartinez-Sansigre et al . (2006) . Such steep slopes are associated 
with evolved (i.e. lobe-dominated) radio sources (IMartmez-Sansigre et al.ll2006l ). while flat 
spectra are generally core-dominated sources. Given our radio resolution, and the redshifts 
of these sources, we are only sensitive to extended structures of > 100 kpc. Below we discuss 
the two sources with observed radio jets. 



3.3. Extended radio emission 



In § 13. 2\ we discussed that our spectral indices suggest evolved radio sources. MIPS15880 
shows a transition FRI/FRII structure with a dia meter of ~ 200 kpc (s ee Figure [2]), which 



is near the mean size for double- lobed structures (Ide Vries et al.l 120061 ) . Given that such 
moderate luminosity radio sources have hot spot a dvance speeds in the range ~ 0.01 -0.1c, 
this size implies a radio source age of ~ 10 7 yrs (e.g. IWillott et al.ll2002l ). 



The bottom panel of Figure [2] (MIPS16122) shows the characteristic tadpole shape 
indicative of a one-sided jet, suggesting that the jet axis is close to the line-of-sight (<< 45°). 
It spans ~ 9Q/sin(i) kpc, where % is the inclination angle. Higher resolution radio data is, 
however, required to confirm the morphology of the radio source. We s ee no counterpart to 
the extended emission in our HST NICMOS image (jPasyra et al.l 120071 ). 



3 Defined as a=log(S 6 ioMHz/Si. 4 GHz)/log(610MHz/1.4GHz). 
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3.4. Radio vs. mid-IR emission 

Figure [3]puts our sources in context with related populations. A number of assumptions 
were necessary to put the various disparate samples on the same plane (see caption Figure [3]). 
These assumptions are crude, but acceptable given the dynamic range of this figure. 

Does our classification agree with the more common 'radio- loud' definition, R = f5GHz/f4400A > 10? 
Given the significant obscuration of our sources this definition is not immediately applica- 
ble. We translate f44 00A to f!4 /im luminosities based on the average quasar spectrum 



([Richards et al.ll2006bl ). and assume a = -0.7 in the radio, which leads to an equivalent cri- 
terion: = Li_4GHz/ Lufim > (dotted line in Figure [3]). Four of our sources remain 
'radio-loud' in this definition as well, while three are borderline, although possibly still 
radio-loud given the rough derivation above. In addition, the screen extinction correction we 
apply to Li^ m is likely to overestimate the mid-IR luminosities. With no strong dichotomy 
between our 'radio-loud' and 'radio-quiet' sources, are we not subject to rather arbitrary 
definitions? Obviously, we see extended emission among the radio- loud sources and not the 
rest of the sample. Only four of the non-radio-loud sources have detections at 610 MHz, 
implying a generally flatter radio spectrum, than for the radio-loud sources. Lastly, the 
surprising difference between high-r and low-r sources (see §3.1). These support the view 
of our radio-loud sources as a distinct population. 



3.5. Black hole masses 



The total IR luminosities (~ 3-4Li 4Atm ; Sajina et al. 2007 in prep.) are likely ~Lboi 
for such highly obscured sources. Based on these, and assuming Eddington accretion, we 
estimate black hole masses of ~2-6x 1O 8 M . This is supported by th e empirical rela- 
tion between radio luminosity, Mbh and accretion rate (ILacy et al.l 1200 ll ). which leads to 
(Mbh) ~ 3.0 x 10 8 M again assuming Eddington accretion. 

Our sources do not represent a complete census of massive black holes at 2~2, and 
are likely to experience further growth between then and z~0. With this in mind, can we 
reconcile such massive black holes at z ~ 2 with the space density of massive black holes 
today? The densi ty of our z~2 heavily obscured radio-loud sources is ~3 x 10~ 7 Mpc -3 . 
Lauer et al.l (120071 ) estimate the space density of black holes with Mbh > 3 x 10 s M (typical 
of our sources) is ~ 3 x 10 -4 Mpc~ 3 locally. This generous gap between the space densities 
implies we are not overproducing black holes at z ~ 2. 
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Discussion 



Our study for the first time directly demonstrates the existence of deep silicate absorp- 
tion (r9.7~l.-6) radio-loud sources at z~2. This is in s trong contrast w ith z<l radio-loud 
populations where none have r 9 . 7 > 1 fjHaas et al.ll2005l ; lOgle et al.ll2006l ). The high obscu- 
rations of our sources is largely a selection bias, as our sample was selected to be mid-IR 
bright and red (§2). However, high radio fluxes are not part of our selection, and hence the 
40% radio-loud fraction of our high-r 2~2 sources is a real effect. This implies a direct 
or indirect coupling between the radio source and surrounding ISM. IWeedman et al.l (120061 ) 
present a z ~ 2 sample similar to ours, but including a radio selection. A quick analysis of 
their sample yields two sources that are both radio-loud, and have r g 7 > 1, further supporting 
to our results. 

The much greater Li4 Mm /-^i.4GHz ratios of our sources relative to local radio-loud sources 
imply we are observing these sources in the brief window between the birth of the radio source 
and before feedback e ffects have shed the dusty envelope and halted the starburst and/or 
black hole accretion. IWillott et al.l (120021 ) argue that the IR-brightness is anti-correlated 
with the age of the radio source consistent with a radio source-driven feedback scenario. 
Given the age of MIPS15880 (>10 7 yrs), its radio jets are not necessarily powerful enough 
to accomplish this. But this may be a function of radio luminosity. MRC1138-262, a. z = 2.2 
radio galaxy more than two orders of magnitude more radio-luminous than our sources, 
does show powerful outflo ws, which are probably driven by lateral shocks from the radio 
jets ( INesvadba et al.ll2006l ). A mid-IR study of a sample of z~2 radio galaxies, covering a 
range of radio luminosities, would therefore provide very useful constraints on radio jet/ISM 
interactions at the epoch when radio galaxies are in the process of shedding their dusty 
envelopes. 

The above does not explain why our sources should have deeper silicate absorption than 
local radio-loud populations. The answer hinges on the unknown origin of the mid-IR ob- 
scuration - s.a. dusty torus or host galaxy. The deep silicate features of our sources imply 
optically thick obscuring media with high filling factors (i.e. likely puffed-up geometry, but 
not A n given the radio jets) as they have to e ffectively 'hide' the hot inner regions of the 
AGN ( IPier fc Kroliklll992l ; iLevenson et al.ll2007f). AGN tori are believed to be clumpy, which 
precludes high fil ling factor s ( Elitzur fc Shlosmanll2006l ). consistent with their obs erved weak 



silica te features (jShi et al.l 120061 ). Higher r 9 .7 can be due to a cold dust screen (lOgle et al. 



20061 ). such as expected from a dusty host galaxy. High- z radio galaxies /rad i o-loud quasars 



have been detected in the sub- mm indicating cold dust (IWillott et al.ll2002l ; iReuland et al. 



20031 ). We already know that the number density o f deeply obscured sources is about 10 
times higher at z ~ 2 than at z ~ (jSajina et al.ll2007l ). Lastly, MIPS16122 is consistent with 
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a jet axis close to the line-of-sight (< <45°). The dusty torus o f an AGN is believed to be 
coupled with the accretion disk (e.g. lElitzur fc Shlosmanll2006l ). implying a nearly face-on 
torus for MIPS 16 122, which is unlikely given its deep silicate absorption feature. On the 
other han d, a mis-match betw een the orientation of the host galaxy and AGN is frequently 
observed (ISchmitt et al.ll2002l ). Thus, although the origin of the mid-IR obscuration remains 
an open question, host galaxy obscuration is a good candidate. 
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Table 1. Summary of radio- loud sample properties 



MIPSID 


z 


T9.7 


l0g(LpAH,7.7) 


log(Li4 Mm ) a 


5l.4GHz 


^lOMHz 


a 


l0g(Li.4GHz) 










L & 


mJy 


mJy 




WHz" 1 


8327 


2.48 


2.4 


10.62 ±0.21 


12.40 ±0.04 


1.40 ±0.06 


3.44 ±0.14 


-1.1 


25.89 ±0.03 b 


15880 c 


1.68 


3.9 


10.60 ±0.09 


12.51 ±0.01 


0.60 ±0.03 


1.35 ±0.06 


-1.0 


25.02 ±0.14 


16059 d 


2.43 


2.7 


10.76 ±0.09 


12.46 ±0.02 


0.57 ±0.03 






25.26 ±0.04 


16122 


1.97 


1.8 


10.41 ±0.21 


12.13±0.09 


2.82 ±0.14 


4.34±0.12 


-0.5 


25.66 ±0.003 


22204 


2.08 


1.6 


11.01 ±0.10 


12.63 ±0.01 


1.89 ±0.08 


4.67±0.10 


-1.1 


25.82 ±0.004 


22558 


3.2 


>5.6 


10.96 ±0.12 


12.88 ±0.01 


0.62 ±0.03 


1.55±0.12 


-1.1 


25.80 ±0.02 


22277 


1.77 


1.4 


10.44 ±0.19 


12.36 ±0.05 


1.57±0.07 


4.79±0.11 


-1.3 


25.65 ±0.005 



a This is monochromatic vL v . 

b This error accounts for the flux error alone, ~0.2dex error in a leads to an additional ~0.1dex 
uncertainty. 

c The radio fluxes for MIPS15880 are the sum of the three components listed in the respective cata- 
logues, we have checked that there is no significant extended flux missed. 

d MIPS 16059 is undetected at 610 MHz, since it falls in the noisy edge of the map. We assume a = -0.7 
for this source. 
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Fig. 1. — Top: The 1.4 GHz luminosity vs. the 9.7/zm optical depth. For our sam- 
ple, only sources with strong silicate features (rg 7 >l) meet the radio-loud criterion 
{Liaguz > 10 25 WHz _1 ). The pink box highlights the region occupied by our seven radio- 
loud sources. Bottom: Mid-IR continuum luminosity vs. T9.7. After extinction correction, 
Lufxm of the high-r sources are higher than those of the low-r sources, possibly providing 
clues as to why their radio luminosities are higher as well. 
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Fig. 2. — The VLA 1.4 GHz contours (1-5 a) overlaid on the MIPS24 /im image for 
MIPS15880 (top), and MIPS16122 (bottom) respectively. At the redshifts of both sources 
the radio beam (~5") subtends Ri42kpc. None of the extended structures have counterparts 
in the i?-band or IRAC3.6 fim images. On the right hand side are shown the IRS mid-IR 
spectra of both sources. In both cases the redshifts are derived from the silicate feature. 
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Fig. 3. — A comp arison between our sources and other radio-loud or ULIRG populations. 
Ogle et al.l (120061 ) quote L 15 , which is close enough to L 14 for our purposes. The Type 2 
quasars we plot are all at z ~ 0.4 - 0.8 so that the observed 24 /zm flux is close t o the rest-frame 



14 nm . The HzRGs shows are those with both 24 fim and 70 /im detections (jSeymour et al. 



20071 ). we interpolated between these fluxes to derive the rest-rame 14 /xm. Both for our 
sample and the ULIRGs we use extinction corrected L 14Mm . The 1.4 GHz fluxes of the PG 
quasars, Type 2 quasars, and ULIRGs are taken from NED. The 1.4 GHz is derived by 
assuming a = -0.7. The dotted line represents the classic radio-to-optical ratio definition of 
'radio-loud' (see text for details). 



